Thermal Equilibrium with the Wiener Potential: 
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Abstract. - We consider the statistical mechanics of a classical particle in a one- 



dimensional box subjected to a random potential which constitutes a Wiener process 
on the coordinate axis. The distribution of the free energy and all correlation functions 
of the Gibbs states may be calculated exactly as a function of the box length and tem- 
perature. This allows for a detailed test of results obtained by the replica variational 
approximation scheme. We show that this scheme provides a reasonable estimate of the 
averaged free energy. Furthermore our results shed more light on the validity of the con- 
cept of approximate ultrametricity which is a central assumption of the replica variational 
method. 



The replica trick and Parisi's scheme of replica symmetry breaking within mean field 
theories provided a promising step towards the application of powerful methods of field 
theory to systems with quenched random couplings [[]]]. More recently these concepts 
have been applied to a variety of physical systems as an approximation scheme within 
a framework called "replica variational approximation" (RVA). The strategy of the RVA 
is to calculate variational bounds on the moments of the canonical partition sum 
within an ensemble of quenched disordered couplings. For integer n > these bounds 
can be achieved from standard inequalities of equilibrium statistical mechanics. Disorder 
averaged physical quantities are estimated from these bounds by using the u n — > 0" trick 
ln(x) = lim n _ (x n — l)/n. The method has been tested using toy models of manifolds in 
random media 0, |J. In the limit of a manifold embedded in infinite dimensional space the 
RVA with a Gaussian trial distribution becomes exact. For a zero dimensional manifold 
in one dimensional space the test has used exact results on certain linear combinations 
of moments of the manifold's position which have been obtained some time ago f|. In 
a more recent work || the probability distribution of the free energy of a related model 
defined on the half axis has also been obtained. 

The purpose of the present letter is twofold. First, we give a complete solution of 
the zero dimensional toy model. We use the expression "complete" to indicate that 
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the solution does not only provide an explicit analytical expression for the probability 
density of the free energy for arbitrary temperatures and system sizes but also allows 
us to calculate all correlation functions of the Gibbs state (canonical distribution). In 
this way we completely characterize the probability functional of the Gibbs state which 
contains all the statistical information on the thermostatics of the disordered system. 
Second we use our results to test the quality of the RVA. 

We consider the position x of a classical particle in a one-dimensional box < x < L 
subjected to a random potential. On the one hand we will use for the random potential 
the standard Wiener process W(x), that is, the Gaussian Markov process with averages 
(W(x)) = and (W(x)W(y)) = mm(x, y). But mainly we will focus our interest on 
its variant V(x) := W(x) — W(L)/2. Below we will refer to these potentials as the 
W- and ^/-ensemble, respectively. A single particle moving under the influence of the 
potential V(x) may be considered as a continuum model of a kink in an Ising chain with 
random fields uncorrelated at different sites |7J whereas the IV-ensemble is related to the 
asymptotic behaviour of a directed polymer in a random potential in 1 + 1 dimensions 
||. Note that the Gaussian random potential V(x) is equivalently characterized by the 
translationally invariant averages (V(x)) = and (V(x)V(y)) = {L/2 — \x — y\)/2. The 
position x of a particle in contact with a heat bath of fixed inverse temperature (3 is 
distributed according to the Gibbs state £lp t z,(x) := exp{— (3V (x)} / Z ^(L) . In the following 
we will be interested in statistical properties of this state, the corresponding partition sum 
Zp(L) := J L dx exp{— (3V(x)} and the free energy Fp(L) := — f3~ l \nZp(L) within the V- 
ensemble. Our general strategy will be to calculate such properties from corresponding 
properties of the VF-ensemble with partition sum Zp(L) := J L dx exp{— /3W(x)} and free 
energy Fp(L) := -/T 1 In Z P (L). 

Throughout our calculations we find it convenient to use scaled variables for locations, 
potentials and partition sums defined as: s := j3 2 x/8, v(s) := — (j3 / 2)V (8s / (3 2 ) , w(s) : = 
-(p/2)W(8s/p 2 ), z{s) := (p 2 /2)Z f3 (8s/p 2 ), and z{s) := {f3 2 /2)Z p {8s/ (3 2 ). The scaled 
system length will be denoted by / := [3 2 L/8. The transition density p s (w, z\wo, z ) of 
the homogeneous Markov process (w(s), z(s)) may be determined from the Fokker-Planck 
equation 

d s p s (w, z\w , z ) = (d 2 w - 4e 2w d^j p s (w, z\w , z ) (1) 

with initial condition po(w, z\wq, z ) = 5(w — wo) 5(z — zq). This equation is equivalent to 
the Langevin equations dw(s)/ds = \^2^(s) and dz(s)/ds = Ae 2w ^ following immediately 
from the definitions of w(s) and z(s). Here £ denotes standard Gaussian white noise with 
= an d (£( s )£(y)) = — s ')- To obtain the joint probability density of v (s) and 
z(s) from the transition density of w(s) and z(s) we use the relation 



S(z(s) — z) 5(v(s) — v) 
= J dwe™| dzpi_ s (w,z v + \w,e w z^jp s (v + \w,e w z 0,0) . 

Note that according to ([!]) one has 

p s (w, z\w , z ) = p s (w, z - z \w , 0). (3) 

To solve the Fokker-Planck equation ([!]) we perform a Laplace transform with respect 
to i and introduce 

roc 

p s (w, g\w ) := dze~ ze p s (w,z\w ,0) (4) 
Jo 
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which obeys the equation 



[d s -d 2 w + Age 2w ) p s (w, q\w q ) = (5) 

with initial condition p (w,g\wo) = 8(w — wo). Thus we see that p becomes a Green 
function of a Bessel-type differential equation. It may be given in terms of an eigenfunction 
expansion 



p s (w, g\w ) = -2 

7T 



/"°° dz/ z/sinh(vrz/) e" 1 ^ K iv {2^~ge w ) K iv {2^e w °) , (6) 

•/ 



where K iv (y) = / °° dt cos(tz/) exp{— ycosht} denotes the modified Bessel function of 
third kind with index iu, see || Chap. 9]. As to the validity of we remark that 
e~ u2s KiJ2^ge w ) obeys (§) for all v which can be checked with the help of the differential 
equation || Eq. 9.6.1] for the modified Bessel functions. Furthermore, the right-hand side 



of (g) reproduces the demanded initial condition which can be inferred from the pair [K| 
p. 173] of reciprocal formulas for the Kontorovich-Lebedev transformation. 

In the sequel we will summarize some of our results on the statistics of the free energy 
and the Gibbs state. All of these results can be obtained from the fundamental solution 



(j^). A detailed exposition of the calculations and the results will be given in [H . 

For the discussion of the statistical properties of the free energies it is convenient to 
introduce their scaled variants f(l) := — Z -1 / 2 ln(z(Z)) and /(/) := — Z -1 / 2 ln(z(Z)) from 
which the original quantities can be recovered according to 

and analogously for Fp(L). As z(l) = e~ w ^z(l), the averaged free energy is the same for 
both ensembles, that is, (/(/)) = (/(Z)). Fluctuations of the free energy, however, will 
differ between the two ensembles. For example, one has ((/(Z)) 2 ) = ((/(Z)) 2 ) — 2 as shown 

3- 
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The density of the free energy /(/) turns out to be 

W(l) - /)) = ^ expj^ + M} K (2erf) 

-— -2e fVI cosh(t) I . 
In Fig. Ill a plot of the centered density of f(l) is given for different values of Z. An 



expression analogous to @ for the density of /(/) will be given in [|Tl[ . 

Of course, the densities allow for a calculation of the averaged free energies and their 
fluctuations. Here we only give the asymptotic expansion of the averaged free energy for 
Z — ► oo 



(/(/)) = </ (i)> = __^ + ^ + o( r>), (9 ) 

where 7 = 0.5772 . . . denotes the Euler-Mascheroni constant |], Eq. 6.1.3]. 

In order to get along with the calculation of the averaged free energy within the RVA 
using the standard Gaussian trial distributions we extend the system to the whole real 
axis by softening the walls. More precisely, we use the potential v(s) + (1 — 2s/Z) 2M , s 
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real, M = 1, 2, 3, . . . , instead of v(s) for the calculations and eventually perform the limit 
M —> oo. For this setting we have re-done the RVA along the lines of ||. The details 
will be given in [[TI]]. The results of the RVA to the averaged free energy (f(l)) ~ f ByA {l) 
are as follows. For I < \^3tt the saddle-point equations have a replica symmetric solution 
only, resulting in 

/RV A (0 ;= /M(i) ;= - (1 + InfSvr/ 2 ^ 2 )) , for I < v^, 

V 71 " 2V/ v v 77 ( 10 ) 

where < V[ < 1 solves + V\ — 1- For / > v^3~7r a replica symmetry broken solution 
exists, too. Since this solution is probably the stable one, we set 

3(1 + V2) 

-= i o t v ^ — mi 071 i i , lur i v 071 . 

mi 



/ RVA C) - / RSB « - - if, ( 3 + ^ - b K» • for ' > ^ 



In Fig. 2 the exact and approximate averaged free energy are compared. One finds 
that the RVA constitutes a lower bound. This has also been observed in comparison with 
simulations [[3], To our knowledge, this phenomenon has not yet been satisfactorily 
explained. 

Now we turn to the calculation and discussion of the statistical properties of the Gibbs 
states. As can be seen from its definition, Vlp ^x) is invariant under a change from the V- 
to the VF-ensemble. Therefore we may calculate all the correlation functions of the Gibbs 
state directly within the IF-ensemble. The J-point correlation function of the scaled 
Gibbs state u>i(s) := (8//9 2 ) ^/3,sz//3 2 (8s/ /3 2 ) can be expressed due to (||) in terms of the 



Fig. 1. - Centered density of the scaled free energy f{l) of the ^/-ensemble for different system 
sizes I = 1/10,1,10,100. 
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Lapace transformed transition density p as 

( II ^( s i) ) = 77 — 7T7 / dgg- 7 " 1 / &w x ■ • ■ dw J+1 Y[p Sj+1 _ s .(w j+1 , g\wj) e 2W] , 
\j=i I l J - L ) ] Jo J ~°° 3=0 (12) 

where =: s < s± < s 2 < ■ ■ ■ sj < sj + \ := I and Wo : = is assumed. With the help of 
(|^) this expression may be simplified further. Eventually we end up with a representa- 
tion of the J-point correlation function in terms of a ( J + l)-fold integral of elementary 



transcendental functions, see HIT 



Here we specialize to the most interesting situation of a large system I ^> 1 and large 
separations between the intermediate points Sj. As to be expected, the leading behaviour 



J 1 J 



limZ J ( YluiQcTj)) = — . _(T i)' O<0i,...,<7j < 1, 

\j=l I 7iy<7i(l - <T X ) j=2 (13) 



l— >oo 



of the correlation function reflects the fact that the system possesses a unique ground 
state with probability 1. For < o\ < a 2 < ■ ■ ■ < oj < 1 the next-to-leading correction 
turns out to be 



hm (Uunilaj)) = i (1 '! )(1 ", )/2 11^ ~ n J: ) **. 

W / irMl-aj) {J-ly- j=2 



(14) 



Fig. 2. - The main plot shows the RVA Fa (L) to the averaged free energy (Fp(L)) for L = 1 
as a function of temperature 1//?. The vertical dashed line indicates the temperature, below 
which replica symmetry breaking sets in. For 1/(3 near both curves are dominated by the 
logarithmic contribution in (0). Therefore the inlets show the scaled variants / RS (Z), / RSB (Z) 
and (/(£)}> respectively, as a function of 1/1. Although the leading correction of the RVA to its 
ground-state energy for I — ► oo is of correct order Z -1 / 2 , it is clearly seen that it has the wrong 
sign. 



5 



A short distance expansion given in [j i\ reveals that the correlation functions are analytic 
for small Sj — Sj-\ and finite I. 

These results allow for a discussion of the validity of the RVA on a deeper level then 
does (H) because we may now check the complete structure of the probability functional of 
Gibbs states. Note that the average (o>(s)) will always be a Gaussian within the replica 



variational ansatz used in || whereas it is actually (tc^Js(1 — s))" 1 for large systems. 
One might argue that this is far away from a Gaussian. However, the Gaussian form is 
inevitably encoded in the ansatz and thus we do not consider this a serious defect. Much 
more interesting is a comparison between the statistics of typical distances A = s—s'. This 
has been obtained in || within the RVA. The result implies that the probability density 
-P(A) should decrease like a Gaussian for large A. The exact result can be obtained from 
(lu(s)lu(s')) and is given by P( A) ~ (l/Vl67r)|A|~ 3 / 2 for large but finite I and A. Thus the 
RVA severely underestimates the range of correlations. However, the most crucial feature 
of the RVA, which is the hierarchical structure of the ansatz, seems to be in accordance 
with a real property of the system, at least in a qualitative sense as we will now explain. 
The geometry of typical distances may be studied by considering higher order correlations. 
For example, consider 3 points < si < s 2 < S3 < I with corresponding distances Ai 2 , A 13 
and A23. If we assume / and all distances to be large but finite, their joint probability 
density is given by P(A 12 ,A 13 ,A 23 ) ~ (l/327r)(A 12 A 23 )- 3 / 2 5(A 12 + A 23 - A 13 ). Thus for 
fixed Si and s 3 the intermediate location s 2 will most probably be near to one of these 
points. If we regard these three points as a "one-dimensional triangle" we conclude that 
isoscele triangles with a shorter third side are statistically preferred. In this sense, there is 
a weak, statistical form of ultrametricity present in the system. An ultrametric geometry 
also appears as a consequence of the hierarchical ansatz for the replica symmetry breaking 
scheme jy, [| which thus captures an important feature of the probability functional of 
Gibbs states. 
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